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Multiple Scattering: Dispersion, Temperature 
Dependence, and Annular Pistons 



Kimball A. Milton, Jef Wagner, Prachi Parashar, Ines Cavero-Pelaez, Iver Brevik, 
and Simen A. Ellingsen 



Abstract We review various applications of the multiple scattering approach to the 
calculation of Casimir forces between separate bodies, including dispersion, wedge 
geometries, annular pistons, and temperature dependence. Exact results are obtained 
in many cases. 



1^ . 1 Quantum Vacuum Energy 

Quantum vacuum energies, or Casimir energies, are important at all energy scales, 

<^ , from subnuclear to cosmological. Applications are starting to appear in nanotech- 

^S| ' nology. Furthermore it is most likely that the source of dark energy that makes up 

VO i some 70% of the energy budget of the universe is quantum vacuum fluctuations. In 

^T ' particular, the 7-year WMAP data is completely consistent with the existence of a 

cosmological constant |[T], 
l> 
O; w=- = -1.10±0.14(68%CL), (1) 
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which is precisely what would be expected if dark energy arose from this source 
llJl . Finally, zero-point fluctuations may be the most fundamental aspect of quantum 
field theory. 



2 Multiple- Scattering Formulation 

The multiple scattering formulation is easiest stated for a scalar field, which is rather 
'easily' generalized to electromagnetism. For example, see Ref. |i3J. Vacuum energy 
is given by the famous trace-log formula. 



E = -TrlnG^ -Tr In GG.', 



where in terms of the background potential V, 
Now we define the T -matrix. 



(2) 



(3) 



(4) 



and if the potential has two disjoint parts, V ~V\+ V2, it is easy to derive the inter- 
action between the two bodies (potentials): 



£■12 = -2Ti"ln(l-GoriGor2) 
= -^Trln(l-yiGiy2G2), 



(5a) 
(5b) 



where G,- = (1 +G0V,) Go, / = 1,2, and likewise IJ refers to V,-. 



3 Quantum Vacuum Energy — Dispersion 



Perhaps not surprisingly in retrospect, we find that the usual dispersive form of the 
electromagnetic energy |4| 



-^M)£S 



'-^EHrHH\r) 



(6) 



must be used, which, quantum mechanically, corresponds to the vacuum energy 

form 

de 
2etrr + to^trr , (7) 



i [ , , , [ do) 



d(0 
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in terms of the Green's dyadic F. This result follows directly from the trace-log 
formula for the vacuum energy 



i f dco 
S = - — Trlnr, 

2J 2n 

and is equivalent to the variational statement Q 



2j 2n 



(8) 



(9) 



From the energy, precisely because the dispersive derivative terms are present, 
we recover the Lifshitz formula for the energy per area between parallel dielectric 
slabs, with permittivity £12, separated by a medium of permittivity £3 of thickness 
a, 

J = ^ J^d^£dkk [In {I -nEr'TEe'^"'") +ln(l -rxMrWe"'"'")] , dO) 



with fc, = \/k]_ + C,'^£i, C, = —i(0 being the imaginary frequency. The TE reflection 
coefficients are given by 



K^-Ki I K-3 - K-2 



K^ + Ki 



K-3 + K-2 



(11) 



while the TM coefficients are obtained from these by the substitution K^ — > K^ = 
Ka/Ea- For further details of this calculation, see Ref. L6J. 



4 Noncontact Gears 

The program of calculating the quantum vacuum lateral force between corrugated 
surfaces and gears has been under active development. The electromagnetic situa- 
tion of currugated dielectric slabs is illustrated in Fig.[T] For details see Ref. [TJ. 



Fig. 1 Parallel dielectric slabs 
with sinusoidal corrugations. 
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In the conductor limit (£,■ — > oo) and for the case of sinusoidal corrugations de- 
scribed by hi (y) — hi sin[A;o(3' + 3'o)] and h2{y) — h2sin[k()y] the lateral force can be 
evaluated to be in first order in hi /a and /12/fl 



where 



(lA) 



ito) = 



fB^^: 



15 



Tt'* J- 



2koa sin(^o3'o) 


p{0) 


In 

a 


a 


dt / sds . , 
Jo sinhi 


sinh 


s+ 


ri (.2 + 4-^2)21 

[2 ^ 8^2^^ J 



(12) 



(13) 



where ,52 ~ s-^ +t^ and s\_ = s^ + {t + to)^. The first term in Eq. ( fT3] l corresponds 
to the Dirichlet scalar case lH), which here corresponds to the E mode (referred to 
in Ref. JU as the TM mode). We note that A'e4li0) = 1. See Fig. |2]for the plot of 
^£^00(^0*2) versus koa. We observe that only in the proximity force approximation 
limit koa = is the electromagnetic contribution twice that of the Dirichlet case, 
and in general the electromagnetic case is less than twice that of the Dirichlet case. 
This result can be shown to coincide with the expression found in Emig et al. ||9l 



versus k^a. The dotted curve 
represents 2 times the Dirich- 
let case. 




apart from an overall factor of 2, which presumably is a transcription error. The 
double integral representation in Eq. ( fT3l ) is more useful for numerical evaluation 
than the single-integral form given in Ref. |9| because of the oscillatory nature of 
the function sinx/x in the latter. Generalization of these results are forthcoming. 



5 Wedge as generalization of cylinder 



In a series of papers, we have considered variations on the wedge geometry, such as 
a wedge defined by perfectly reflecting walls, intersected with a concentric circular 
cylinder, the arc being either a perfect reflector itself, or the boundary between two 
dielectric-diamagnetic regions. Most interesting is the case when the wedge itself is 
constructed as the interface between two such media. See Fig. [3] In order to have a 
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tractable situation, we have considered the diaphanous or isorefractive condition 



£1^1 = £2^2, 



(14) 



that is, the speed of light is the same in the two media. (If that is not done for the 
wedge, the differential equations are no longer separable.) See Refs lfT0l[mri2i for 
more detail. 



Fig. 3 Wedge geometries, 
(a) The perfectly conducting 
wedge geometry, (b) The ge- 
ometry of a wedge intercut by 
a perfectly conducting cylin- 
drical arc. (c) Wedge with 
magnetodielectric arc. (d) Di- 
aphanous wedge in a perfectly 
conducting cylindrical shell. 





£i,Vii 
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m~^ 


X 



Consider now case (d). Using multiple scattering, or the Kontorovich-Lebedev 
transformation, we obtain the following implicit formula for the eigenvalues for 
the order v of the contributing cylindrical partial waves, D{v,Co) = 0, where (r = 
reflection coefficient on wedge) 



D{V,(0) = (l_e2;i,V)2_^2(g/v(2;i-a)_g,-va)2 

= -4e2'^'^[sin2(v7r) - r^sm^{v{7r - a))], 



(15) 



which are selected by the "argument principle," which is just the Cauchy theorem 
applied to the contour 7 shown in Fig. HI 

In this way, we find the energy per length given by <^ = „ ' i ^{p), P — 1^1 CL as 



shown in Fig.|5] Note that only for perfect reflectors does the energy diverge as the 
opening angle approaches zero. 



6 Annular Piston — Semitransparent Plates 



The wedge geometry may be generalized by considering two semitransparent plates 
in a Dirichlet annulus, as shown in Fig. |6] We use multiple scattering in the angular 
coordinates, and an eigenvalue condition in the radial coordinates; this problem is 
equally well solvable with radial Green's functions, but this approach may be more 
generalizable. This section is based on Ref. II13I . 



Fig. 4 Contour of integration 
7 for the arguinent principle. 
Sliown also are singularities 
of the integrand along the real 
and imaginary v axes. 
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Fig. 5 The function e{p) plot- 
ted as a function of opening 
angle a. 
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Fig. 6 Two semitransparent 
plates in an annulus. 




The Green's function {#(r,r') will satisfy the equation 

[-V2- a)2 + y(r)] ^(r,r') = 5(r-r'), 



(16) 



while '^^ ' has V (r) = 0. For the cylindrical geometry of an annulus, the boundary 
conditions are ^ = at p = a and p = b, where a and b are the inner and outer 
radii, respectively. We take the potential to be V (r) —v{9)/p^. The corresponding 
Green's function is 



dk 



(r,r';a))=/ ^e'*(")^7;^(p;«,^)«„(p';a),%^(0,0'), (17) 
/-oo zn ^ 
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in terms of the separation constant 77 . The normaUzed radial eigenfunctions appear- 
ing here are 



dp dp 



Rr^{p;co,k)^ri^Rr^{p;(0,k), (18) 



with the boundary conditions Rri {a;(0,k) = Rfj {b; (0,k) = 0. The reduced Green's 
function satisfies 



-^2W+m 



8^(9,9') = d{e-e'), 



(19) 



with periodic boundary conditions. 

To obtain the radial functions, we need the solution of the modified Bessel dif- 
ferential equation, of imaginary order, which is zero for p ~ a for all values of 77 
and K. An obvious solution is 



Rr^{p;K) ^ Kir^iKa)Iir^{Kp) -Iir^{Ka)Kir^{Kp) ^ R^t^{P,k) 



(20) 



where 



Iv^-{Iv+I- 



(21) 



The eigenvalues are given by the zeros of Z)(tj) = Rrj{b;K). We don't need the 
explicit eigenfunctions here. 



6.1 Reduced Green's Function 



The free angular reduced Green's function is given by 



gf{9,9') = ^(-^xr^i)\9-9'\ + '^^co,h^\9-9'\]. (22) 



For a single potential v{9)—X5{9 — a) for 9,9' <E [a,27Z + a], the reduced Green's 
function is 



„ 1 / ,„ „/, 277cosh77;rcosh77|0 - 0'| 

gni0,9') ^ — -sinhT] 0-0' +-!——! ^-^i- ^ 

'^ ^ 277 V 277smh77;r + Acosh77;r 

cosh 77 (2;r + 2a - - 0') - cosh277;rcosh 7] 1 - 0'| 

— A } 7 ^ I . (23) 

[277 sinh 77;r + A cosh 777rJ2sinh777r 
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6.2 Two Semitransparent Planes 

Now we look at the interaction energy between two semitransparent planes, as il- 
lustrated in Fig. |6] Since it is nontrivial to work out the Green's function for two 
potentials, it is easiest to use the multiple-scattering formalism dSbl) 

E^l r ^Trln(l -^('Vi$^P)y^)_ (24) 

The subscripts on the Vs represent the potentials yi(r) = Xi5{0)/p^, and V2(r) = 
^25(0 — a)/p^. The Green's functions with superscript (/) represent the interac- 
tion with only a single potential Vj. From this we obtain a simplified form of the 
interaction energy: 

^ = ^ i" ^dKj^in (l - tr^^'vi^^Vz) , (25) 

where ^i are given in Eq. ( |23] l. Then 

tr ^'V ^^V == AiA2Cosh^T7(7r-«) 

'' '' (277sinh777r + AicoshT77r)(277sinh777r + A2Cosh777i;) 

Using the argument principle to determine the angular eigenvalues, we get the 
following expression for the energy for an annular Casimir piston, 

1 /"OO /" I) 

S' = —^. / KdK / drj— In [Kij^ {Ka)Iin {Kb) - Iij^{Ka)Kir^ {xb)] 
7 
rAnfl Mhcosh^rjJTZ -a)/ cosher] 7Z \ 
\ (2TjtanhTj;r + Ai) (2TjtanhTj;r + A2)/ 

The contour of integration for the argument principle is again given in Fig. |4] 

This formula can actually be used to evaluate the energy of interaction between 
the two planes of the piston, by distorting the rj contour to lines making angles of 
±k/4 with respect to the real axis. The results are shown in Fig. [7] In Fig. |7] we 
define d = ^^ sin ^ , and the plateaus seen in the second figure may be understood 
from the proximity force approximation. 



'g'pFA ^ 1 b^ 

S\\ 16 



^0 + ^) ' (28) 



in comparison to the interaction between infinite parallel plates. 
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Fig. 7 Energy/length for an annular piston as function of angle (top), and compared to the en- 
ergy/length for a rectangular piston (bottom). 



7 Applications of Multiple Scattering 

As an illustration of practical calculations using the multiple scattering machinery, 
we illustrate in Fig.|8]a semi-infinite array of periodic potentials, such as a array of 
dielectric slabs, for which the exact Casimir-Polder force with an atom to the left 
may be calculated lfT4l . 



V(x) 



Fig. 8 A semi-infinite an'ay 
of periodic potentials. The 
exact CP force between an 
atom and this array may be 
calculated. 



d A 2A 3A 4A 5A 6A 
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7.1 Casimir-Polder Force 

Consider an atom, of polarizability a{(i)), a distance Z to the left of the array. The 
Casimir-Polder energy is 






E = -l dC I j:^^a{iQtig{Z,Z), (29) 



where apart from an irrelevant constant the trace of the Green's function is 

trg(Z,Z) ^ — [-C2^TE^(^2^2fe2)^TM]^-2^|Z|^ (3O) 

Here the reflection coefficients are those for the entire array (a is the distance be- 
tween the potential slabs), 



^ = 1- 
2R 



„2Ka I n2 t2 



R^-T^-J {e^^" -R^-T^y- 4R^T^ 



(31) 



If the potentials consist of dielectric slabs, with dielectric constant e and thick- 
ness d, the TE reflection and transmission coefficients for a single slab are (k"' = 

V^C^TF) 

^^^ = ^^-7^ ., ,, . , (32a) 



( i+k'/k \ 2K'd _ (h:JSUl\ 



TE _ 4{K'/K)e''''' 

The TM reflection and transmission coefficients are obtained by replacing, except 
in the exponents, k' -^ k' /e. (Multilayer potentials have been discussed extensively 
in the past, see, for example, Refs. lIBllTtl lTTllTSll .) 

For example, in the static limit, where we disregard the frequency dependence of 
the polarizability, 

E = -^^F{a/Z,d/Z). (33) 

This is compared with the single slab result in Fig. |9] It is interesting to consider the 
Z — > 00 limit, which is shown in Fig. [TO] When a/d -^ we recover the bulk limit. 
Such results apparently will have applications to experiment rather soon 1191 . 



8 Exact Temperature Results 

The scalar Casimir energy between two weak nonoverlapping potentials Vi (r) and 
V2(r) at temperature T is I.20J 
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Fig. 9 Casimir-Polder energy 
between a semi-infinite an'ay 
of dielectric slabs with e = 
2, compared to the energy 
(lower curve) if only one 
slab were present. Here we 
have assumed that the spacing 
between the slabs and the 
widths of the slabs are equal. 



Fig. 10 Casimir-Polder en- 
ergy for large distances from 
the array, as a function of the 
ratio a/d, where a is the dis- 
tance between the dielectric 
slabs in the array, and d is the 
thickness of each slab. Here 
e = 2. 




Ej 



32n^ 



{dr){dr')V,{r}V2{r') 



r'|2 



(34) 



8.1 Exact Proximity Force Approximation 



From Eq. ( |34] | we find that the energy between a semitransparent plane and an arbi- 
trarily curved nonintersecting semitransparent surface is for weak coupling 



Ej — — 



X1X2T 



dS I dx , 

2nTz{S) X 



(35) 



where the area integral is over the curved surface. Here z{S) is the distance between 
the plates at a given point on the surface S. Equation (l35T l is precisely what one 
means by the proximity force approximation: 



(36) 



EpFA = J dS^ll{z{S)), 
as noted by Decca et al. I.21J . See also Ref. 
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8.2 Interaction Between Semitransparent Spheres 



We can, for weak scalar coupling, compute the energy between two spheres of radius 
a and b, whose centers are separated by a distance R: 

X1X2 ab( l-la-bf/R^ , , .. , , 

~f{2nT{R-a + b))-f{2KT{R + a-b))y 
where f{y) for y < TT is given by the power series. 



m = L 



„fi2«(2«-l)(2«)! 
which is obtained from the differential equation 



■/", 



(37) 



(38) 



d 



y^ifiy) = '^°*3' - -, /(o) = /'(o) = 

ay^ y 



(39) 



Results for the energy obtained by solving this differential equation are shown in 
Fig.[TT] For further details see Ref. 



Fig. 11 Comparison between 
the general and high temper- 
ature forms of the energy, as 
a function of RT. Energies 
are shown for a = b = R/4. 
The high temperature re- 
sult is linear in T. Also 
shown is the power series 
expansion truncated at 200 
terms, which diverges in this 
case atRr = 1/3. Plotted is 
e = -16;rR£'/(AiA2«2). 




8.3 Mean distances between spheres 



Encountered in the above calculation are mean powers of distances between spheres 
as defined by 



dndn'\r -r'\P ^ {4-7:)^RPPp{a,b), 



(40) 



Multiple Scattering 



13 



for spheres, of radii a and b, respectively, separated by a center-to-center distance R. 
Here a = a/R and b = b/R, and Pp(a, b) can in general be represented by the infinite 

series 

2 r(2n-p-\) 



Here the homogeneous polynomials Q„ are 

Go = 1, 



-Q„{d,b). 



Qi = 3a^ + lOcP-P + 3b'', 



S21.4 



23 = 4fl*' + 28a^/7^ + 28fl^/7* + 4Z7 



(41) 



(42a) 
(42b) 
(42c) 
(42d) 



Here in general. 



e« = xL 



1 A /2n + 2 



There is also a recursion relation. 



\2m + 1 



R-P d 



'2(ji-m)r2m 



Vi(«,^) = Y^^7;'+%(«,^) 



(43) 



(44) 



since 2„ is homogeneous in R of degree — 2n. 

For integer p > —2, Pp is a polynomial of degree 2\p/2], and we can immedi- 
ately find 



PpiaM 



1 



1 



[il+a + b)P+^ 



4ab{p + 2){p + 3) 

+ {l-a-bf+^-{l-a + b)P+^-{l+a-b)"+^], (45) 



Although this was derived for integer p it actually holds for all values of p. 

For example, when p is a negative integer, we have the explicit forms, which are 
obtained from Eq. ( |45l l by taking the appropriate limit: 



P_i = 1, Newton's theorem, 



1 

P-2 = -^ 

Aab 



, l-{a + bf „, (1+i 
In ^ J-+a\n- 



)2-fl2 



P_3 = -TT^ln 

Pa 



\-{a~bf 

\-{a + bf 
Aab 1 ~ (fl - ^)2 ' 
1 



;i-/7)2-fl2 



1 



[\-{d + bf][\-{c 



'b\n 



(l+a)2-P 

(l-fl)2-^2 



(46a) 
(46b) 

(46c) 
(46d) 



and further expressions, which can be obtained by use of Eq. (l44l i. may be readily 
verified. 
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9 Conclusions 

The multiple scattering formalism can be used to find numerical results effectively 
in many situations, as we have seen in this outline. Weak coupling results are exact 
and often given in closed form. The method can also be used to extract not only 
interaction energies but self energies, as described in Ref. ll23l . 
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